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Abstract 

We calculate the spin correlation function and the magnetic longitu- 
dinal and transverse susceptibilities of a two-dimensional antiferromagnet 
doped with a small concentration of holes, in the t-J model. We find that 
the motion of holes generates spin fluctuations which add to the quan- 
tum fluctuations, the spin correlations decaying with the inverse of the 
spin distance, while increasing with doping as the critical hole concen- 
tration, where the long-range order disappears, is approached. Moreover, 
the longitudinal susceptibility becomes finite in the presence of doping, 
due to the strong damping effects induced by the hole motion, while the 
transverse susceptibility is renormalized by softening effects. Both the 
longitudinal and the transverse susceptibilities increase with doping, the 
former more significantly than the latter. Our results imply that doping 



destroys the long-range order while local antiferromagnetic spin correla- 
tions persist. This is consistent with experiments on the doped copper 
oxide superconductors. 

PACS: 71.27.+a, 74.25.Ha, 75.10.Nr, 75.30.Cr 
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Since their discovery/ the copper oxide high-Tc superconductors have shown 
unusual magnetic characteristics, along with the unconventional transport properties.^ 
The undoped materials, e.g., La2Cu04, are antiferromagnetic (AF) insulators, 
and doping, e.g., in La2-5Sr5Cu04, introduces holes, which are the charge car- 
riers, in the spin lattice of the copper oxide planes. The Cu02 planes are 
described by a spin- 1/2 Heisenberg antiferromagnet on a square lattice with 
moving holes that strongly interact with the spin array. A remarkable feature 
of the copper oxides is the strong dependence of their magnetic properties on 
the hole concentration 6. In previous work^~^ we studied the effects of doping 
on various magnetic properties, and showed that the motion of holes generates 
significant softening and damping of the spin excitations, leading, in particular, 
to the disappearance of the long-range AF order at a small hole concentration, 
due to the decay of spin- waves. We found that the staggered magnetization 
vanishes at a hole concentration well below the one for which the spin-wave 
velocity vanishes, or even the one for which all spin-waves become overdamped. 
This suggests that although the long-range order has disappeared, strong AF 
correlations persist, which allow the spin- wave excitations to exist. This is in 
agreement with experiments in the copper oxides, which show that, although 
the long-range order disappears, AF correlations persist up to fairly high dop- 
ing, into the superconducting state. It is therefore of interest to study the 
spin correlations in these materials, because of their unusual behavior and their 
possible connection to high-Tc superconductivity. 

In this work we use the t—J model to calculate the spin correlation function 
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of a two-dimensional antiferromagnet as a function of the hole concentration, 
which allows to investigate the local spin fluctuations, and also calculate the 
longitudinal and transverse magnetic susceptibilities, which reflect the global 
response of the system, accounting for the total spin fluctuations. We consider 
zero temperature and the low doping regime where the long-range AF order 
still exists. It is shown that the motion of holes generates spin fluctuations 
which add to the quantum fluctuations of the system, and increase with the hole 
concentration. Moreover, we flnd that the longitudinal spin susceptibility, which 
is zero in a pure Heisenberg antiferromagnet at zero temperature, becomes flnite 
in the presence of doping, increasing significantly with the hole concentration, 
more pronouncedly than the corresponding transverse spin susceptibility. 
We describe the copper oxide planes with the t-J model, 

Ht-j = -t {cl^Cja + h.c.) + J (Si-S-j-^n^nA, (1) 

<i,j>,cr <i,j> 

where Sj = |cJ^(Ta/3Ci/3 is the electronic spin operator, a are the Pauli matrices, 
m = riii + riii and riia- = cl^Ci^- To enforce no double occupancy of sites, 
we use the slave-fermion Schwinger boson representation^" for the electronic 
operators Cia = fjha, where the slave-fermion operator fj creates a hole and 
the boson operator bia- accounts for the spin, subject to the local constraint 
fifi + ^It^iT + ^ii^il — ^-S*- For the undoped system, the model (1) describes a 
spin- 1/2 Heisenberg antiferromagnet, exhibiting long-range Neel order at zero 
temperature. The Neel state is represented by a condensate of Bose fields bi^ = 
^/2S and bji — V2S, respectively, in the up and down sub-lattices, and the 
bosons bi = bi^ and bj = are then spin-excitation operators on the Neel 
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background. After a Bogoliubov-Valatin transformation on the boson Fourier 
transform, bk = Ukf3k + Vkl3^_^, where Wk = [((1 — 7k)~^^^ + l) /2] and 
Vk = -sgni^jk) [((1 - 7k)~^^^ ~ l) /2]^^^, with 7k = {coskx + cosky)/2, we 
arrive at the effective Hamiltonian 

^ = - ^ E /q4-k [^(q' -k)/3-k + V{q - k, k)/?t] + ^a;£/3^/3k, (2) 



q,k 



where, y(q, k) = zt{'-fqUk + Tq+ki'k) represents the interaction between holes 
and spin- waves resulting from the motion of holes with emission and absorption 
of spin-waves, wJJ = (zJ/2)(l — 7k)^^^ is the dispersion for spin-waves in the 
undoped antiferromagnet, and z is the lattice coordination number (z = 4), A'' 
the number of sites in each sub-lattice. The sums are performed in the first 
Brillouin zone of an antiferromagnet on a square lattice. 

The magnetic properties are calculated in terms of the spin-wave Green's 
functions 

D-+{k,t- t') = -i< T(3k{t)pl{t') > 
D+-{k,t- t') = -i< T0l^{t)0-k{t') > 
D—{k,t - t') = -i <T/3k{t)l3^k{t') > 
D++{k,t- t') = -i< Tf3l^{t)pl{t') >, 
where < • • • > represents the average over the ground state. The spin-wave 

Green's functions verify the Dyson equations 

D^^ (k, to) = Z?r (k, ^) + E (k, a;)n"^ (k, (k, to) , 

with /X, r;,a, 7 = =b. The free Green's functions are Z)~+(k, a;) = l/(a; — + 
iv), D+-{k,u;) = l/{-uj-co^ + irj), (r, ^ 0+), D--{k,uj) = D++{k,u;) = 



0. n"'>'(k, w) are the self-energies generated by the interaction between holes 
and spin-waves, which we calculate in the self-consistent Born approximation 
(SCBA). This corresponds to considering only "bubble" diagrams with dressed 
hole propagators, describing the decay of spin- waves into "particle- hole" pairs. 
The spin-wave self-energies take the form* 

(k, w) = 1 ^ (k, q) [F(q, -k; a;) + F(q - k, k; -w)] , (3) 
q 

with [/— (k,q) = [/++(k,q) = y(q,-k)y(q-k,k), f/+-(k,q) = y(q-k,k)2, 
C/-+(k,q) =l^(q,-k)2 and 



/'+00 pO 

y(q,-k;a;) = / duj' du" 

Jo J-oo 



CO + CO" — CO' + IT] 

The SCBA provides a spectral function for the holes, p(q, w), which is com- 
posed of a coherent quasiparticle peak with weight tto — (J/t)^/^ and dispersion 
£q — £mm + (q— qi)^/2"^> with effective mass m ~ 1/J, the Fermi surface for the 
holes consisting of pockets, of radius qp = Vtt5, located at q^ = (±^,±|) in 
the Brillouin zone, and an incoherent continuum taking the approximate form 
he{\u\ - zJ/2)0{2zt + zJ/2 - \uj\), with /i ~ (1 - ao)/2zt. We calculated the 
self-energies to lowest order in the hole concentration i5. 
The spin correlation function is defined as 

= ^ E ((S,- • S,-+r) - (S,-) • (S,+r)) , (4) 



where the sum runs over all lattice sites. Writing the spin operators, Sj, SJ = 

{S^ + Sj)/2, Sj = {S^ — Sj)/2i, in terms of the electron operators, one has, 
•^1 ~ (4t'^-?T ~ ^J.'^JJ-)/^' '^t ~ '^h^jl' '^i — 4i^-'f ' "wliich, using the Schwinger 



boson representation and the boson condensation associated with the Neel state, 
leads to 5? = (1 - 5){l/2 - b]bj), Sf = (1 - S)bj, Sr = (1 - S)b], for the up 
sub-lattice, and = -(1 - 6){l/2 - bibj), S+ = (1 - d)b], S' = (1 - 6)bj, for 
the down sub- lattice, having done the approximation /j/j = 1 — 6. In terms of 
the spin-excitation boson operators one has 



C(r) = (1-6) 



E 



4N 



b]b^+^)+2[b]bjb]^,b,+, 
[(b]b,+r) + (b,bl 



jes(i) 



+2 (6t6,6t^^6,+r) - 2 (6t6,) {b\^^b,+r)) 



(5) 



After Fourier transform and the Bogoliubov-Valatin transformation, we make 
the mean-field decoupUng {ABCD) « {AB) {CD) + {AC) {ED) + {AD) {BC). 
This allows to express the correlation function (5) in terms of the spin-wave 
Green's functions in the form 



C(r) = {-lY+y{l-5)^\^^Y.^ul+vl)cos{).-v) 



(6) 



p-\-oo 

/ -^{21mD+-{k,uj)+Au]^vwlmD--{k,u))) 
Jo 27r 



— V 



+ 00 



duji r°° (hj2 



m /V Jo 27r 7o Stt 

Kl,K2 



cos((ki - k2) • r) 



X [vl, - {ul, + vlJlmD+ {ki,uji) -2uk,VkJmD (ki,wi)] 
X [ul, - (wL + vl^)lmD+-{-k2,0J2) - 2uu^vu^ImD—{k2,W2)] } , 

where r = {x,y). The pre-factor of (—1) arises when the correlation is between 
sites on different sub-lattices. 
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To lowest order in the hole concentration, S, we obtain for the correlation 
function ©, the expression 



C{r) = (1-5)^ [Coir) + Cs{r)] 



(7) 



where 



Coir) 



^ Xl(^k + ^^k) cos(k • r) 



(8) 



^ (Y1 cos(k • r) j vl cos(k • r) j 



iV2 

is the correlation function for a pure Heisenberg antiferromagnet, accounting 
for the quantum fluctuations, and 



Csir) 



1 + ^ II("k + ^^k) cos(k • r) cos(k • r) 



(9) 



WkWk 



Ren++(k,t^°) „ 



+ « + vl) 



+ 00 



+ 00 



dwlmn++(k,w) 



Imn-+(k,t^) 



contains the effect of doping on the spin correlations associated to the hole mo- 
tion; the pre-factor (1 — S)^ corresponds to spin dilution, being negligible in 
the low doping regime considered. In Fig.l we present the correlation function 
C(r), Eq.(|ZI), calculated for two different directions, x — y and y = 0, in the 
case oi 6 — 0.02, and the pure case, S = 0.0. One sees that the spin correlations 
are independent of the spatial direction, a result that is verified at any doping. 
In Fig. 2 we plot the correlation function as a function of the spin distance r, 
for various hole concentrations. We observe that C(r) increases with doping, 
and decays, at large distances, as 1/r (Fig. 2 inset), both in the pure and in 



the doped cases. One can describe the behavior of C(r) at large distances as 
C(r) = A{5)/r, where ^(,5) = -\- B5", with = 1/V2n and a = 0.42. 
A{d) contains the doping dependence, which is represented in Fig. 3. The dom- 
inant contribution to Cs{r), at large r, comes from the imaginary part of the 
spin- wave self-energies, which depend on the hole concentration essentially as 
V^.^ In Fig. 4 we compare the increase of C(r) with the hole concentration S 
at fixed small r (Fig. 4a) and large r (Fig. 4b). The decay of C(r) with 1/r was 
expected for the undoped case since in a two dimensional antiferromagnet at 
zero temperature the correlation length is infinite. ^'^^ One also expected C(r) 
to increase with doping since the motion of holes generates spin fluctuations 
which eventually lead to the destruction of the long-range AF order at a fi- 
nite critical concentration 6c- In previous work,^ we found that the staggered 
magnetization vanishes at a small critical concentration (e.g. 6c — 0.07 for 
t/J = 3), while the long- wavelength spin excitations remain well defined up to 
a higher hole concentration (i5* ~ 0.17 also for t/J = 3). Here, we find that 
the doping does not qualitatively change the behavior of C(r) with r, as com- 
pared to the pure case, which reflects the robustness of the local AF order in 
the doped material. Spin correlations in the copper oxides were studied before, 
both experimentally^'^"^ and theoretically, ^^"^o ^ higher doping regime 

where the long-range AF order has already disappeared. In this regime the spin 
correlations decrease with increasing doping, as the system moves away from 
the critical hole concentration. 

In the presence of long-range AF order one distinguishes a longitudinal and 
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a transverse susceptibility. The longitudinal spin susceptibility is defined as 



XII =X||(k = 0,u; = 0), (10) 



where the dynamical susceptibility is given by 

r'+oo 



X|| 



/•+00 

(k,w)=i / dte''^* {[S'{k,t),S'{-]si,0)]) . 
Jo 



In terms of the spin-wave Green's functions one has 



XII ^ ^ 



k^o^N ^ / 2^ [2wki^^kiUki-kWki-k - Wki^ki-k - ^'ki^'ki-k] 
x£)+-(ki,a;)£)-+(ki -k, -w) 



which to lowest order in the hole concentration gives 



1 /■+°°rfa; Imn+-(k,a;) 



The transverse spin susceptibility is defined by 

X± = X±(k = 0,a; = 0), (12) 

where 

X±{k,oj)=i / dte^-'{[S^(k,t),S^i-k,0)]}. 
Jo 

In terms of the spin-wave Green's functions the transverse spin susceptibility is 
expressed as^ 

/I _ \ 

= - lim [Re£>+- (k, 0) + Re£>++ (k, 0)] , 

k^O V 1 + 7k / 

which, to lowest order in the hole concentration 5, is given by 

2 



X_L = lim 



k^o zJ{l + 7k) 



1 - ^2)1/2 [ReH"*" (k. 0) + Ren++ (k, 0)] 

(13) 
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We found that X-L takes the simple form 

x± = z^xl 

where Xj_ = l/C^zJ) is the transverse spin susceptibiUty for a pure Heisenberg 
antiferromagnet, and Z^^ = 1 + 4(5a^ (j)^ is a renormalization factor. 

Comparing Eqs. (11) and (13) one sees that the motion of holes influences 
the longitudinal and transverse susceptibilities in different ways, the former is 
produced by the imaginary part of the self-energy, while the latter is renormal- 
ized by the real part of the self-energies. In a pure Heisenberg antiferromagnet 
the longitudinal susceptibility is zero. However, with doping x\\ acquires a finite 
value due to the decay of spin- waves into "particle-hole" pairs, generated by 
hole motion. The renormalization of X-L reflects a softening of the spin cou- 
pling induced by the hole motion. In Figs. 5 and 6 we plot the longitudinal, 
Eq.(ll), and the transverse, Eq.(13), susceptibilities as a function of the hole 
concentration for t/J = 3, 4, in the doping range where the long-range AF order 
exists, in the approach considered.^ We find that both susceptibilities increase 
with doping, though the longitudinal one is far more sensitive to the hole con- 
centration than the transverse one. The transverse susceptibility reflects the 
stiffness of the antiferromagnetic lattice. In contrast, the longitudinal suscep- 
tibility is set by the strong damping effects, which are also responsible for the 
disappearance of the long-range AF order at low doping.^ When the long-range 
order is broken, the susceptibility of the system should be essentially given by 
X. ~ |x|| + |X-L) with the longitudinal susceptibility providing an important 
contribution. Also, in the ceramic samples whose crystal axis are randomized 

11 



the susceptibility x is given by an average of the susceptibilities for the three 
directions. An increase of the spin susceptibility with doping has in fact been 
observed experimentally. ■^^"^^ 

In summary, we studied the effects of hole motion on the spin correlation 
function and the magnetic longitudinal and transverse susceptibilities of a two- 
dimensional antiferromagnet doped with a small concentration of holes. We 
found that the spin fluctuations increase with doping, the spin correlations 
decaying with the inverse of the spin distance, which indicates that the local AF 
correlations remain quite robust. Furthermore, we show that the longitudinal 
magnetic susceptibility acquires a finite value in the presence of doping, due to 
the strong damping effects generated by the hole motion, while the transverse 
magnetic susceptibility is renormalized. Both susceptibilities show a significant 
increase with doping, which is however more pronounced in the longitudinal one. 
Our results imply that doping destroys the long-range AF order while local spin 
correlations persist. This is consistent with experimental observations in the 
copper oxide high-T^ superconductors. 

*Electronic address: fdias@cii.fc.ul.pt 
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Figure Captions: 

FIG. 1. Correlation function C(r) vs spin distance r, at the hole concentration 
S = 0.02, for the directions x = y (open circles) and y = (diamonds) on the 
square lattice, with t/J = 3. Inset: C(r) vs r, in the pure antiferromagnet 
{S = 0.0). 

FIG. 2. Correlation function C(r) vs spin distance r, x = y, for various hole 
concentrations 6, with t/J = 3. Inset: C(r) vs 1/r for, large r. 

FIG. 3. Spin correlation amplitude A{S), at large spin distances, vs hole con- 
centration S, with t/J = 3. 

FIG. 4. Correlation function C(r) vs hole concentration ^, with t/J = 3, for 
fixed r, X = y. a: in the range of small r. b: in the range of large r. 

FIG. 5. Longitudinal susceptibility x\i ^ function of doping 6, for t/J = 3 
(open circles) and t/J = A (diamonds). 

FIG. 6. Transverse susceptibility X-L as a function of doping 5, for t/J = 3 
(open circles) and t/J = 4: (diamonds). 
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